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ABSTRACT 

Using infrared, radio, and 7-ray data, we investigate the propagation characteristics of cosmic-ray 
(CR) electrons and nuclei in the 30 Doradus (30 Dor) star-forming region in the Large Magellanic Cloud 
(LMC) using a phenomenological model based on the radio-far-infrared correlation within galaxies. 
Employing a correlation analysis, we derive an average propagation length of ~ 100 — 140 pc for 
~ 3 GeV CR electrons resident in 30 Dor from consideration of the radio and infrared data. Assuming 
that the observed 7-ray emission towards 30 Dor is associated with the star-forming region, and 
applying the same methodology to the infrared and 7-ray data, we estimate a ~ 20 GeV propagation 
length of 200 — 320 pc for the CR nuclei. This is approximately twice as large as for ~ 3 GeV CR 
electrons, corresponding to a spatial diffusion coefficient that is ~ 4 times higher, scaling as (R/GV) S 
with 6 s=a 0.7 — 0.8 depending on the smearing kernel used in the correlation analysis. This value is 
in agreement with the results found by extending the correlation analysis to include ~ 70 GeV CR 
nuclei traced by the 3—10 GeV 7-ray data (6 sw 0.66 ± 0.23). Using the mean age of the stellar 
populations in 30 Dor and the results from our correlation analysis, we estimate a diffusion coefficient 
Dr « 0.9 — 1.0 x 10 27 (i?/GV) 0,7 cm 2 s _1 . We compare the values of the CR electron propagation 
length and surface brightness for 30 Dor and the LMC as a whole with those of entire disk galaxies. 
We find that the trend of decreasing average CR propagation distance with increasing disk-averaged 
star formation activity holds for the LMC, and extends down to single star-forming regions, at least 
for the case of 30 Dor. 

Subject headings: cosmic rays - galaxies: individual (LMC) - 7 rays: galaxies - Hn regions - infrared: 
galaxies - radio continuum: galaxies - stars: formation 



1. INTRODUCTION 

Cosmic rays (CRs) are a dynamically important com- 
ponent of the interstellar medium (ISM) in galaxies. Yet 
their role in shaping galaxy evolution is currently un- 
clear due to the difficulties in chara cterizing their origi n 
and propagation in the ISM (e.g., iStrong et al.l 120071 ). 
The energy density of CRs is comparable to that of mag- 
netic fields, as well as the radiation fields and turbu- 
lent motions of the interstellar gas in galaxies. Together, 
CRs and large-scale magnetic fields comprise a relativis- 
tic plasma whose interactions with the interstellar gas 
shape the overal l chem i stry and heating of the ISM (e.g., 
IBoulares fc Coxl H99CJ iFerrierd l2lJ0l 1Coxl 12001 . and 
may also play a signific ant role regulating star-formation 
processes (e.g..lSocrates et~al . 2008; Papadop oulosll20iol ; 
iPapadopoulos et al.ll2011f ). 

The diffuse emissions from the radio to high-energy 7 
rays (> 100 MeV), produced by various interactions be- 
tween CRs and the interstellar gas, radiation and, mag- 
netic fields, are the best way to characterize the physics 
of CRs throughout most of the Milky Way, and other 
galaxies. For example, the CR lepton component can be 
probed by radio, X-ray, and 7-ray observations. The ra- 
dio emission arising from the injection of CR electrons 
from supernova remnants trace these particles as they 
propagate through large-scale magnetic fields and lose 



energy from synchrotron radiation. Similarly, the same 
CR electrons inverse Compton (IC) scatter off the in- 
terstellar radiation field (ISRF) and cosmic microwave 
background, yielding o bservable e missions from X-ray to 
7-ray energies (Po rter et al. 2008). High-energy 7 rays 
are particularly useful because this energy range gives 
access to the dominant hadronic component in CRs via 
the observation of 7 rays from the decay of neutral pions 
produced by inelast ic collisions between CR nuclei and 
the interstellar gas (jPollack &: Faziolll963h . 

Externally viewed galaxies have the advantage of mit- 
igating line-of-sight confusion that hampers interpreta- 
tion of the diffuse emissions of the Milky Way. While 
such observations can be used to determine the present 
distribution of CRs, understanding their propagation 
history remains difficult without knowledge of the ini- 
tial distribution of CR sources. Using the tight, em- 
pirical correlation between the far-infrared (FIR) and 
(predominantl y) non-thermal radio continuum emiss ion 
from galaxies (|de Jong et al.|[l985t iHelou et al.lll985ft . a 
number of studies have attempted to characterize the 
propagation of CR e l ectrons in external galaxie s (e.g., 
iMarsh fc HeToul fliM iMurphv et al.ll2006al 12001 . The 
underlying physics relating the radio and FIR emission 
from galaxies lies in the process of massive star forma- 
tion: young massive stars are the dominant sources of 
dust heating, and end their lives as supernovae whose 
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Figure 1. The point-source subtracted 1.4 GHz (left) and 24 /an (right) maps surrounding 30 Dor with residual 1 — 3 GeV 7-ray contours 
starting at 0.5 X 10 6 counts sr — 1 (roughly ~1 count per 0.1° X 0.1° pixel) with linear intervals of 0.326 X 10 6 counts sr _1 overlaid. The 
1.4 GHz map has been corrected for free-free emission. The region shown is centered at RA 05 h 39 m 07 s , Dec— 69°22'02" (J2000) and has 
a radius of 1° (si 875 pc). The residua ls between the radio/7-ray and infrared map s were calculated over this region ( H2.40 . as well the 
minimum energy magnetic field f i|2,2.1H , and the radiation field energy density ( 92.41 . In addition, we show the locations of two Crab-like 
pulsars, PSR J0540-6919 (square) and PSR J0537-6910 (circle), along with a known background X-ray source RX J0536.9-6913 ('X') 
that are coincident with the line-of-sight towards 30 Dor on the 1.4 GHz and 24 /im maps. 



remnants presumably accelerate and inject CR elec- 
trons into the ISM where they produce diffuse syn- 
chrotron emission. It was hypothesized that the radio 
image of a galaxy should resemble a smoothed version 
of its infrared image, because the mean free path of 
dust-heating photons (~ 100 pc) is significantly shorter 
than the typical diffusion length of CR electrons (~ 
1 — 2 kpc) (|Bicav fc He lou 1990). This phenomenology 
is su pported by stud i es wi thin nearby galaxies at kpc 
(e.g., Marsh fc Heloul IT998) and few hundred pc (e.g., 
iMurphy et al.l l2006al ) scales, as well as o n scales of > 
50 pc in the Large Magellanic Cloud (LMC . lHugh.es et al.1 
2006). Investigations on sub- kpc scales within a sam- 
ple of nearby spirals have even shown a dependence of 
the typical propagation length of CR electrons on star- 
formation activity arising from the predominant youth 
of CR electron populations in galaxies wit h enhanced 
disk- averaged star formation rates (|Murphy et al. 2006b, 
[2001) . 

Recently, using 11 mon ths of data, the L MC was de- 
tected by the Fermi- LAT ()Abdo et alj2010bl) at high sig- 
nificance with the resolved 7-ray emission showing very 
little correlation with the large-scale distribution of the 
gas column density. If CRs freely diffuse in the ISM of 
the LMC, as they appear to in the Milky Way, the 7- 
ray emission should be correlated with the distribution 
of gas, which is pr edominantly neutral hydrog en and he- 
lium in the LMC (jStavelev-Smith et alj |2003). Instead, 
the observed 7-ray emission for the LMC is more strongly 
correlated with tracers of massive star-forming regions. 
Employing the multi-frequency diffuse emissions from ra- 
dio to 7-ray energies, we can investigate for the first time 
the propagation of both CR electrons and nuclei associ- 
ated with a star-forming region in an externally resolved 
galaxy, 30 Dor in the LMC. 

2. DATA AND ANALYSIS 

We compiled archival infrared and radio data, together 
with 32 months of 7-ray data from the Fermi-LAT for 
our analysis. We assume a distance to 30 Dor of 50 kpc 



throughout this analysis 1 . 

2.1. Spitzer Infrared Data 

Spitzer imaging of the LMC was carried out as part 
of the Surveying the Agent s of a Galaxy's Evolution 
(SAGE: IMeixner et all 120061 ) legacy program. The ob- 
servations covered an area of ~ 7° x 7° over the LMC 
in all IRAC (3.6, 4.5, 5.8, and 8 /im) and MIPS (24, 70, 
and 160 /im) bandpasses. The individual observations 
were calibrated, combined, and mo saicked by the SAGE 
team, using a custom pipeline (see IMeixner et al.l 120061 
for further details). Because we are only interested in 
diffuse emission, we use point-source subtracted versions 
of these images prepared by the SAGE team (see right 
panel of Figure [TJ . For the present study, we only make 
use of the MIPS 24, 70, and 160 /im data, which have 
native resolutions of « 5'.'7, 17", and 38", respectively. 
The photometric uncertainties are conservatively taken 
to be 5, 10, and 15% at 24, 70, and 160 /im, respectively 
(see the MIPS Instrument Handbook 2 ). 

The majority of the emission at 24 /im arises from dust 
heated in the vicinity of massive stars, which are the 
progenitors of core-collapse supernovae. Consequently, 
the 24 /xm data are used as a proxy for the CR source 
distribution in our phcnomcnological modeling. We use 
the full infrared data (i.e., 24, 70, and 160 /jm imaging) 
to estimate the radiation field energy density over the 
LMC and 30 Dor (see ^4l). 



2.2. Parkes+ATCA Radio Data 

We use combined 1.4 GHz single-dish (Parkes 
64 m) and interfero metric (ATCA) data presented in 
iHughes et al.1 (|2006l ). where a much more detailed de- 
scription of the data preparation can be found. We em- 
ploy a version of the radio map that has point sources 

1 Taken from the NASA Extragalactic Database (NED; http:// 
nedwww.ipac.caltech.edu), where there are currently 275 refer- 
ences for individual distance measurements to the LMC 

2 http : //irsa . ipac . caltech . edu/data/SPITZER/docs/mips/ 
mipsinstrumenthandbook/ 
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Figur e 2. The sam e region on the sky as shown in Figure [T] but now displaying the Hi column density map (Stavclcy-Smith ct al. 
2003; Kim ct al. 2003). In the left panel, residual 1 — 3 GeV 7-ray contours are overlaid on the Hi column density map with the same 
intervals as shown in Figure [l] In the right panel, residual 3 — 10 GeV 7-ray contours are overlaid on the Hi column density map starting 
at 0.2 X 10 6 counts sr - * with linear intervals of 0.178 X 10 — 6 counts sr _1 . Also shown are the two Crab-like pulsars, PSR J0540— 6919 
(square) and PSR J0537— 6910 (circle), along with a known background X-ray source RX J0536.9— 6913 ('X') that are coincident with the 
line-of-sight towards 30 Dor. 



removed (see left panel of Figure [l} from the map using a 
median filter technique (see Hu ghes et al. 2006) because 
more than 90 % of such sources are expected to be back- 
ground AGN (|Marx et al.lll997fh The angular resolution 
of the final combined mosaic is 40", and sets the resolu- 
tion for our analysis between the radio and infrared data. 
At our assumed distance for the LMC, this projects to a 
linear scale of s=s 10 pc. 

Because we are interested in the non-thermal com- 
ponent of the radio data, we subtract the contribution 
from thermal (free-free) emission using a scaled version 
of the 2 4 /itm image following Equation 2 of lMurphv et al.l 
(2006b). This relation is based on the empirical cor- 
relation found between the 24 //m and extinction- 
corrected Paa lum inosities from star-f orming regions 
within NGC 5 194 bvlCalzetti et all (|2005|) Although not 
unive rsal (e.g. JPerez-Gonzalez et al.ll2006l : [Calzetti et aD 
2007), it has proved to be a good first-order estimate of 
the thermal radio emission in nearby galaxies when com- 
pared to therm al fractions d erived from multifrequency 
radio data (e.g.. lMurphv et al.ll2008l) . 

To check the reliability of our free-free emission esti- 
mate, we compare with the thermal fraction derived from 
standard radio spectral decompositions usin g single-dish 
radio data from the literature: 1.4 G Hz (jKlein et al.l 
fl989h. 2.3 GHz (IMountfort et alJ ll98l. and 2.45 GHz 
( Havnes et alJll991l ). There is a known discrepancy be- 
tween the total flux from the 1.4 GHz single-dish map 
a nd the combined Par kes+ATCA map being used here, 
as lHughes et al.l (|2007t l report a global 1.4 GHz flux den- 
sity that i s 1.3 times smaller than the single dish esti- 
mate from lKlein et al.l (|1989f ). This difference is largely 
attributed to an underestimation of the beam width for 
the single-dish data. Because the single-d ish data were 
reduce d and analyzed self-consistently by lHavnes et al.l 
(1991), we assume that each map is affected in a similar 
way. 

Using these multifrequency radio data, and assuming 
a constant non-ther mal spectral index of ont ~ 0.7 
(jHavnes et al.l I1991D , we derive a thermal fraction at 
1 GHz that is sa 1.2 times larger than the value of 25% 



found using the 24 /im map over our region of interest 
surrounding 30 Dor shown in Figure [1] This difference is 
consistent with results from a high-resolution r adio con- 
tinuum survey of M 33 (|Tabatabaei et al.ll2007l) , where a 
standard radio spectral decomposition was found to yield 
thermal fractions for 1 1 Hn complexes that were overesti- 
mated by an average factor of 1.5 at 8.3 GHz. However, 
even if we scale the 24 /zm-derived free-free map by a 
factor of 1.2 to match the total thermal fraction mea- 
sured by the radio spectral decomposition in our region 
of interest, our results are not significantly affected (see 
33~2l). 

2.2.1. CR Electron Energy Estimate 

We estimate the typical energies of the CR electrons 
emitting the observed non-thermal 1.4 GHz emission. 
For electrons propagating with a pitch angle 9 in a mag- 
netic field of strength B with isotropically distributed ve- 
locities, such that < sin 2 9 >= | leading to B^ ta 0.825, 
then a CR electron emitting at a critical frequency v c 
will have an energy 
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Integrating the free- free corrected 1.4 GHz flux over the 
regio n shown in Figure Q] (S , ^J GHz = 102.26 Jy; see 
§2.41) . we use the revise d minimum energy calculation 
of iBeck fc Krausg (|2005h . assuming a proton-to-electron 
number density ratio of Kq f=a 100 ± 50, a non-thermal 
spectral index of «nt ~ 0.7 ± 0.1, and a path length of 
!«li 0.5 kpc, to calculate a minimum energy magnetic 
field strength of ~ 11 ± 4 /xG. Using this range of values 
in Equation [TJ our 1.4 GHz maps are most sensitive to 
~ 3 ± 0.5 GeV CR electrons. If the minimum energy 
condition does not hold, and we instead assume an ex- 
treme range in B, such as 3 — 50 fiG, the corresponding 
CR electron energies will range from 6—1.5 GeV. How- 
ever, the (in-)validity of the minimum energy assumption 
is uncertain and often assumed (as we do in this paper). 



MURPHY ET AL. 



Table 1 

Effective 7-ray and Proton Energies for Different 
7-ray Energy Ranges 



Proton Energy 
Index 



7-ray Energy Range 

3 GeV 3-10 GeV 



(GeV) 



(GeV) 



(GeV) 



(GeV) 



2.1 


1.63 


22.4 


5.10 


69.2 


2.2 


1.61 


19.5 


5.04 


60.4 


2.3 


1.60 


17.2 


4.99 


53.6 


2.1 


1.58 


15.4 


4.93 


48.1 


2.5 


1.60 


14.1 


4.87 


43.6 


2.6 


1.55 


12.5 


4.82 


39.8 


2.7 


1.54 


11.5 


4.77 


36.6 


2.8 


1.53 


10.5 


4.72 


33.9 



We briefly discuss the impact of this on th e derived prop- 
erties of the diffusion coefficient in §4.2.11 

2.3. Fermi- L AT Data 
The Fermi-LAT instrument , event reconstruc tion, and 



response are described in lAtwood et al.l (|2009T l. In this 
paper, we use events and instrument response functions 
(IRFs) for the standard low-background "Clean" events 
corresponding to the Pass 7 event selections 3 . To mini- 
mize the contribution from the very bright Earth limb, 
we restrict the event selection and exposure calculation 
to zenith angles < 100°. We selected all Front-converting 
events, for which the point-spread function (PSF) is nar- 
rowest, within a 20° square region of interest (Rol) cen- 
tered on RA 05 h 38 m 42 s , Dec. -69°06'03" (J2000) for 
ss 32 months of sky survey data from 2008 August 4 un- 
til 2011 March 24. Exposure maps and the PSF for the 
pointing history of the observations were generated using 
the standard Fermi-LAT ScienceTools package available 
from the Fermi Science Support Center 4 . The exposure 
of the instrument over the Rol for the data taking period 
used in this analysis is very uniform. 

We constructed foreground-subtracted 7-ray maps 
over the same region using the standard Fermi- 
LAT diffuse emission model 5 together with the ap- 
propriate isotropic background model and 7-ray point 
sources from the Fermi-LAT Secon d Source Catalog 
(|The Fermi-LAT Collaboration! 120111 ). forward folding 
the combination with the exposure and PSF to obtain the 
total counts expected for these foregrounds. These were 
subtracted from the data and the residual emission for 
the 1 — 3 GeV energy interval is shown in Figure Q] over- 
laid on the free-free corrected 1.4 GHz and 24 fim maps. 
The effective 68% containment radius for the Fermi-LAT 
PSF over this energy bin for the event selection used is 
~ 0.50° (« 440 pc), assuming an E~ 2A spectrum for the 
7-ray emission, which is a su fficient approximatio n to the 
spectral shape determined bv lAbdo et al.l (|2010b[ ) for the 
LMC 7-ray emission over this energy range 6 Note that 
the effective full width at half-maximum (FWHM) of the 



3 We use the P7V6.CLEAN event class in this work; see h ttp:// 
fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/ 
Pass7_usage . html 



http : //f ermi . gsf c . nasa . gov/ ssc/data/ analysis/ 

5 http : //f ermi . gsf c . nasa. gov/ssc/data/access/lat/ 
BackgroundModels . html 

6 Weighting with an &JT ' spectrum gives a marginally different 
result, being larger by «3.5%. 



Fermi-LAT PSF over this energy bin is « 0.38° (« 330 
pc), significantly smaller than the 68% containment ra- 
dius quoted above. 

Similarly, we constructed a map for the 3 — 10 GeV 
residual emission. This is shown in Figure [2] where we 
overlay the 1 — 3 a nd 3—10 GeV contours on the H I col- 
umn density data (|Stavelev-Smith et al.ll2003t Kim et aT] 
2003). The effective 68% containment radius for this 
map is ps 0.24° (w 210 pc), with an effective FWHM of 
« 0.19° (« 170 pc), again assuming an E~ 2S spectrum 
for the 7-ray emission over this energy bin. However, 
the 3—10 GeV map has much lower statistics, even near 
30 Dor, contrasting with the 1 — 3 GeV 7-ray map, which 
has > 3 times more events. Therefore, we base the ma- 
jority of our spatial analysis on the 1 — 3 GeV map, and 
only use the higher energy 7-ray data for testing the en- 
ergy dependence of the CR propagation (see i j4.2.2[) . 

2.3.1. CR Nuclei Energy Estimate 

As for electrons and the 1.4 GHz emission, we esti- 
mate the typical energies of the CR nuclei emitting 7 
rays for the energy ranges that we consider in this pa- 
per. For non-AGN dominated star-forming galaxies like 
the Milky Way and the LMC, the assumption that the 
CR nuclei are the predominant relativistic particle pop- 
ulation producing the 7-ra y emission in the F ermi-LAT 
energy range is reasonable (jStrong et al.ll2010l) . 

We define the function 



G(E 7 ,E p ) 




JT7' a dJ P dcr Pp( E i, E p 

dE ^dE " 



E^(E-, 



dE' p P 



dEj 
dJ p da pp {E~ ll E' p ) 



dE p 



dE-, 



(2) 



where E p and E~ are the proton and 7-ray energies, 
/3 = v/c, dJ p /dE' p ex E'~ a is the CR proton flux, 
da pp fdE 1 is the differenti al cross section for 7- r ay pro - 
duction calculated as in Moskalcnk o &: Strong! (1998), 
and -E' nm is the minimum proton energy required to pro- 
duce a photon of energy E^ . The effective proton energy, 
E° , for a given E 7 is given by the root of G(Ej,E p ) so 
that contributions to the 7-ray flux at E 1 below and 
above E p s are equal. The effective proton energy is a 
function of the proton spectral index E p s — E p R (a), 
which we show for different a for the 1 — 3 GeV and 
3 — 10 GeV energy ranges in Table [TJ We also show the 
effective (average) 7-ray energy for the 1 — 3 GeV and 
3 — 10 GeV energy ranges, calculated assuming the same 
spectral index for the 7 rays as for the protons. This as- 
sumption is valid for thin-target interactions by the CR 
nuclei with gas, which applies for the ISM surrounding 
30 Dor and the greater LMC. 

The effective proton energy changes by a factor of 
~ 2 when the power-law index of the proton spectrum 
changes from 2.1 to 2.8, while the effective energy of 7 
rays does not change appreciably. For the 7-ray ener- 
gies, this is not surprising because the effective energy 
is calculated over a relatively narrow energy bin. The 
effective proton energy change is a combination of the 
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fairly steep rise in the inclusive pion production cross 
section as the momentum of the interacting proton in- 
creases, approximately oc E 05 between 10 and 100 GeV 
(see, e.g., Figure 2a from lDermerlll98a ). When the pro- 
ton spectrum is steep (e.g., with an index ~ 2.8), this 
rise is not large enough to compensate for the steeply 
falling number of high-energy protons. When the proton 
spectrum is flat (e.g., with an index ~ 2.1) the rise in 
the inclusive cross section increases the contribution of 
high-energy protons, making the E° s significantly larger. 

2.4. Image Registration and Smearing Analysis 

All images were cropped to a common field-of-view and 
regriddcd to a common pixel scale. To properly compare 
each image at the same resolution, maps were convolved 
with an appropriate PSF. For the radio and infrared (24, 
70, and 160 /im) data, we convolve all images to the 
resolution of the 1.4 GHz image, which has a FWHM of 
40", using a Gaussian beam having a FWHM equal to 
the quadrature difference between the final and original 
FWHM values of each image. 

To match the resolution of the 24 /im data to that of 
the 7-ray maps, we convolve the 24 /im image (at its 
native resolution) with the in-flight Fermi-LAT PSF, de- 
scribed in Section l2~3l Because the FWHM of the 24 /im 
PSF is < 1% of the effective FWHM over the 7-ray en- 
ergy bins used in our analysis, we do not attempt to 
correct for this additional broadening. 

Next, we apply our image smearing analysis (described 
below) to the resolution-matched images to determine 
the extension of the non-thermal radio and 7-ray mor- 
phologies relative to the 24 /im morphol ogy. The pro- 
cedure largel y follows that described in iMurphv et al.l 
(|2006bl l2008f) . where we refer the reader for a more de- 
tailed description. Our underlying assumption is that 
the only difference between these distributions is due to 
the diffusion and energy losses of the CRs because of the 
common origin for the CR leptons and nuclei. We fo- 
cus on the 24 /jm data in the smearing analysis rather 
than the 70 /im data, as was done by IMurphv et al.1 
(e.g.. l2008i ). because it traces warmer dust and has been 
found to be more peake d near star-forming regions than 
the 70 /im morphology (Hclou et al. 2004; IMurphv et al.l 
l2006al) . For instance, typical 24/im knots surrounding 
30 Dor tend to be ~10 pc in size. Thus, the 24 /zm 
maps likely act as a better source function for recently 
accelerated CRs. 

To summarize the procedure, we convolve the entire 
24 /im image of the LMC by a parameterized kernel k(t) 
and compute the residuals between the free- free corrected 
1.4 GHz and smoothed 24 /im maps. The smoothing 
kernel is a function of a two-dimensional position vector 
r, having a magnitude of r = (x 2 + y 2 ) 1 / 2 , where x and 
y are the right ascension and declination offsets on the 
sky, respectively. We investigate both exponential and 
Gaussian kernels because a preference of kernel type may 
suggest different CR transport effects. The exponential 
and Gaussian kernels have the forms of n(r) = €~ T l l and 
ft(r) = e~ r I 1 , respectively, where I is the e-folding scale 
length. Thus, for the Gaussian kernel, a 2 = I 2 /2. 

Exponential kernels, which have broader tails than 
Gaussian kernels of the same scale length, imply energy 
loss and/or escape timescales less than, or comparable to, 



the diffusion timescale. Gaussian kernels, on the other 
hand, suggest a simple random walk scenario. Because 
the LMC lack s a well-definied disk, and has a low inclina- 
tion (i « 31°: lNikolaev et aLll2004f L we assume spherical 
symmetry for 30 Dor. Hence, there is no need to modify 
the position vector by geometric factors. 

We calculate the normalized-squared residuals between 
the radio and the infrared images 



m 






(3) 



where R is the radio image with free-free emission re- 
moved, 1(1) is the infrared image smoothed by a kernel 
with scale length I, Q = HIj/'ERj is used as a normal- 
ization factor (i.e., the global 24 /zm/1.4 GHz ratio), and 
the subscript j indexes each pixel. Only pixels detected 
with a significance of > 3a of the RMS noise in each 
map are used in the calculation of the residuals. The 
minimum in <f> defines the best-fit scale length, which is 
taken to be the typical distance traveled by the CR elec- 
trons. Similarly, we perform the same procedure using 
the 24 /im and 1 — 3 GeV 7-ray maps to determine the 
typical distance traveled by the CR nuclei. 

For this procedure, generically, if there is no broaden- 
ing due to propagation effects, the minimum of Equation 
[3] will be at I — 0, because both the radio and 7-ray data 
have already reduced to the same angular resolution. Be- 
low, we obtain kernel scale lengths indicating that signif- 
icant smoothing of the 24 /im image is required to im- 
prove the match with the radio and 7-ray images. While 
the scale lengths obtained for the radio/24 /tm residual 
analysis are much larger than the angular resolution of 
these maps, for the 7-ray/24 /tm data, the derived scale 
lengths are comparable to the angular resolution after 
convolution with the Fermi-LAT PSF. The technique is 
capable of detecting scale lengths that are a small frac- 
tion of the resolution (i.e., <C AZ = 50 pc, which is the 
scale length step size used in the analysis), and is more 
sensitive for higher signal-to-noise ratio maps. We em- 
phasize that the kernel scale lengths obtained in this case 
are in addition to the smoothing of the 24 /im data to 
match the Fermi-LAT PSF. If there was no additional 
effect from propagation, the derived scale lengths would 
be zero, whereas we detect a meaningfully non-zero scale 
length. 

The uncertainty in <f> is estimated by numerically prop- 
agating the uncertainties in the input images as mea- 
sured by the la RMS noise of each map, with the uncer- 
tainty for the best-fit scale length then estimated as the 
range in scale length corresponding to that from min((/>) 
to min(0) + unc(0) along the residual curve. We note 
that this places a lower limit on the uncertainty for the 
best-fit scale lengths since there is additional uncertainty 
on the 1.4 GHz thermal fraction estimation and fore- 
ground subtraction for the 7-ray maps. The calculation 
of the residuals and photometry was carried out within 
an aperture having a radius of 1° (w 875 pc), encompass- 
ing 30 Dor, centered at 05 h 39 m 07 s , -69°22' 02" (J2000, 
see Figure [T|). 

We use the 24, 70, and 160 /im photometry to calculate 
the total infrared (IR, 8 — 1000 /im) luminosity over this 
region by fitting these data to the spectral energy distri- 
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Figure 3. The residuals between the free -free corrected 1.4 GHz 
and smoothed 24 /im (as defined in i|2,4l) plotted as a function 
of both exponential (asterisks) and Gaussian (diamonds) of kernel 
scale lengths. The shaded region indicates the extrapolated prop- 
agation length for GeV protons based on the compa ct morphol- 
ogy of the 7-ray emission obser ved by the Fermi-h AT IIAbdo et al.l 
2010b). The value reported bv lAbdo efall I2010bft is the a from a 
modeled Gaussian. The values for the Gaussian scale lengths agree 
within errors. 




bution (SED) models of iDale fc Heloul ([2002) and inte- 
grating the best-fit SED between 8 and 1000 ^tm. The in- 
dividual 24, 70, and 160 ^m flux densities are 0.42±0.02, 
4.1 ±0.41, and 8.0 ±1.7 xlO 4 Jy, respectively. In addition 
to photometric uncertainties, the above errors include a 
term for the mean RMS noise values in the convolved 
24, 70, and 160 [iva. maps, which are measured to be 
~0.045, 0.685, and 13.2 MJy sr _1 , respectively. We ob- 
tain an IR luminosity of Ltr = 1.15 ± 0.12 x 10 42 ergs -1 
(3.00 ± 0.31 x 10 8 £ Q ). Taking this value, we estimate 
the corresponding radiation field energy density 



(4) 



where Jbol is the bolometric surface brightness, c is the 
speed of light, and Atr ~ 2.4 kpc 2 is the area over which 
the photometry was measured. All quantities in Equa- 
tion |4] are in cgs units. This calculation is for radiation 
emitted near the surface of a semitransparent body, and 
the parenthetical term provides an empi rically der ived 
correction for non- absorbed UV emission (Bell 2003), re- 
sulting in a value of C/ ra d ~ 2.39 ± 0.17 x 10~ 12 erg cm" 3 
(1.49 ± 0.11 eVcm -3 ) averaged over the volume consid- 
ered. 

3. RESULTS 

We present our results from comparing the morpholo- 
gies of 30 Dor as measured by warm dust emission at 
24 /jm, synchrotron emission at 1.4 GHz, and 7-ray 
emission at 1 — 3 GeV. We compare the morphologies 
both qualitatively, through a visual inspection, as well 
as quantitatively using the image smearing analysis de- 
scribed above. 

3.1. 24 (J-m, 1.4 GHz, and 1 — 3 GeV 7 -ray 
Morphologies 

Figure Q] shows the residual 1 — 3 GeV 7-ray emission 
overplotted as contours on the free-free corrected 1.4 GHz 
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Figure 4. The residuals between t he 1 — 3 GeV 7-ray and 
smoothed 24 (im images (as defined in i|2.4H plotted as a function 
of both exponential (asterisks) and Gaussian (diamonds) of kernel 
scale lengths. The asterisk and diamond attached to the horizontal 
dotted lines indicate the location of the best-fit scale length mea- 
sured for the CR electrons using exponential and Gaussian kernels 
(see Figure [3}. 

radio and 24 /im warm dust emission in the left and right 
panels, respectively. While there is a general correspon- 
dence between the 7-ray, radio, and infrared emission, 
the peak of the 7-ray emission is clearly offset from the 
peak of the radio and infrared emission, which appear 
to be nearly co-spatial. This lack of correspondence be- 
tween the peak of the 7-ray emission with the radio and 
infrared emission leads one to question whether the ob- 
served _^rav_emissionJs_Jndeed associated with 30 Dor. 
In Abdo et all (|2010bl ). it was noted that 7- 
ray emission near 30 Dor may have a non- 
negligible contribution from two Crab-like pulsars, 
PSR J0540-6919 (05 h 40 m 11.2 s , -69°19' 54"; J2000) 
and PSR J0537-6910 (05 h 37 m 47.4 s , -69°10' 20"; 
J2000). In addition, there is a background X-ray source 
(RX J0536.9-6913: 05 h 36 m 57.8 s , -69°13' 26"; J2000) 
that is also coincident with the line-of-sight toward 
30 Dor, which could also contribute to the observed 
emission. The locations of these point sources are shown 
in Figure [1] and Figure [2j and are near, but not at, 
the peak of the 1 — 3 and 3 — 10 GeV 7-ray emission. 
It is also interesting to see that the 1 — 3 and 3 — 10 
7-ray contours appear to peak close to a region of 
high Hi column density in Figure [U It is possible to 
subtract fitted 7-ray point sources at the the locations 
of the pulsars and AGN, resulting in a situation where 
30 Dor does not exhibit 7-ray emission in excess of the 
subtracted foreground model. The implications of such 
a scenario are described in §4.1. However, we note that 
there is no detection of pulsed 7-ray emission so far 
reported from either of the pulsars. Also, no variability 
is detected from the region during our period of analysis. 
If RX J0536.9-6913 is a 7-ray emitting AGN, the lack 
of variability makes it difficult to determine if this 
candidate source contributes significantly in our Rol. 

3.2. Results from the Smearing Analysis 

In Figure [3] we plot the residuals (c/>, see H2.4p between 
the free-free corrected 1.4 GHz and smoothed 24 /zm 
maps as a function of exponential and Gaussian ker- 
nel scale lengths. The residuals between the free- free 
corrected 1.4 GHz and 24 /im maps are decreased by 
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more than a factor of ^2 after smoothing the 24 /im 
map using either exponential or Gaussian kernels. A 
slightly larger improvement is found by using an expo- 
nential kernel, suggesting energy losses and/or escape 
may be important for the CR electrons. The correspond- 
ing best-fit exponential and Gaussian kernel scale lengths 
are lOOtJo and 200^5°° pc, respectively. We note that 
the estimated ~3 GeV CR electron propagation length 
reported here assuming random walk diffusion (i.e., using 
a Gaussian kernel) is consistent within errors to the GeV 
CR pro ton confinement length reported by lAbdo et al.l 
(|2010bD . whose estima t e assum ed Gaussian profiles (see 
Figure [3]). lAbdo et all (|2010bl) report the Gaussian a to 
be a — 170 ± 60 pc, so we multiply this number by \[2 
for proper comparison with the scale length definition for 
our Gaussian kernel. 

As discussed in ^2.2[ we test how our results are af- 
fected by assuming that the thermal fraction of the 
1.4 GHz radio map is being underestimated by the 24 /im 
maps relative to what is derived using single-dish ra- 
dio data at 1.4, 2.3, and 2.45 GHz. By scaling the 
24 /im-derived free- free maps by a factor of 1.2 to match 
the total thermal fraction measured by the radio spec- 
tral decomposition, and repeating the smearing analysis, 
we find that the best-fit exponential scale length is still 
100 pc, while the best-fit Gaussian scale length is slightly 
increased to 250 pc. A larger improvement is again found 
using an exponential kernel relative to a Gaussian ker- 
nel. The results therefore do not appear to be signifi- 
cantly affected by increasing the thermal fraction esti- 
mate to match that from the radio spectral decomposi- 
tion method. 

Similarly, in Figure^ we plot the residuals between the 
1 — 3 GeV 7-ray and smoothed 24 /im maps as a function 
of exponential and Gaussian kernel scale lengths. The 
residuals between the 1 — 3 GeV and 24 /im maps are 
decreased by a factor of ^2 after smoothing the 24 /mi 
map using either exponential or Gaussian kernels. While 
a slightly larger improvement is found by using exponen- 
tial kernels in the comparison between the 1.4 GHz and 
24 /im maps, we do not find a preference in kernel type 
for the 1 — 3 GeV 7-ray map. The corresponding best-fit 
exponential and Gaussian kernel scale lengths are 2001-[q 
and 450^50 pc, respectively. These best-fit scale lengths 
are significantly larger (i.e., a factor of ~2) than what was 
measured for the 1.4 GHz maps. Because the 1.4 GHz 
and 1 — 3 GeV 7-ray maps probe different energy CR par- 
ticle populations, the corresponding differences in their 
best-fit scale lengths may arise from different diffusion 
speeds. This scenario is discussed in ^4.21 

4. DISCUSSION 

Using a phenomcnological image smearing model, we 
have estimated the typical propagation lengths of ~ 
3 GeV CR electrons and ~ 20 GeV CR nuclei by compar- 
ing the spatial distributions of smoothed 24 /im maps im- 
ages with (free- free corrected) 1.4 GHz and (foreground 
subtracted) 1 — 3 GeV 7-ray maps, respectively. From 
the analysis, there is a difference between the typical dis- 
tances travelled by ~ 3 GeV CR electrons and ~ 20 GeV 
protons (assuming a proton energy index of 2.1), where 
the 20 GeV CR protons are found to travel ~2 times fur- 
ther, on average, independent of kernel type. However, 



it is additionally worth including a brief discussion of the 
physical implications for the case where 30 Dor may not 
be emitting in 7 rays. 

4.1. The Case for a Dark 30 Dor in 7 rays 

As stated in §3.1, there are both pulsars and a back- 
ground X-ray emitting source that may contribute to the 
observed 7-ray emission towards 30 Dor. For this case, 
the situation has 30 Dor emitting at 1.4 GHz, therefore 
containing diffusing CR electrons, but no observational 
signature for the CR nuclei. Explaining this would re- 
quire the CR nuclei to escape the system without inter- 
acting with the interstellar gas. 

Taking the above estimate for the bolometric lu- 
minosity from 30 Dor (i.e., the IR luminosity given 
in JJ2.4I corrected by the parenthetical term in Equa- 
tion [4] for non-absorbed UV emission) together with the 
up dated star formation rate (SFR) calibrations given 
in IMurphv et al.l (J2011D . we estimate a correspond- 
ing SFR of ~0.15 Moyr^ 1 . This value is consistent 
with others in the literature (e.g.. IHughes"c t al. 2007; 
IHarris fc Zaritskvl I2009t ILawton et al.l 120101 ) . Assum- 
ing this value and the sensitivity of the Fermi-LAT 
data used in our analysis, the absence of 7-ray emis- 
sion from 30 Dor clearly disagrees with the empirical 
correlations between SFR or the product of supernova 
rate and the gas mass with 7-ray luminosi ty found for 
the local group and nearby st arbursts (e.g. JAbdo et al.1 
I2010al iLenain fe Walterll2011h . Furthermore, theoretical 
expectations of such scaling relations and their implica- 
tions for the contribution by normal star-forming galax- 
ies to the diffuse 7-ray extragalactic back ground (e.g., 
iPavlidou fc Fields! l200ll iFields et al.l [2010h would then 
be questionable. However, because there is currently no 
strong evidence linking the 7-ray emission from 30 Dor 
to the known point sources in the field, we discuss the 
physical interpretation of our analysis with the assump- 
tion that it arises from 30 Dor. 

4.2. Cosmic Ray Propagation from 30 Dor 

To interpret our findings on the kernel scale lengths, 
we must consider the dependence of gas density on the 
synchrotron and 7-ray emission. The diffuse 7 rays in 
the energy range that we have used (1 — 3 GeV) are 
predominantly from CR nuclei interacting with the in- 
terstellar gas, and therefore S~ t oc n p ?iisM where n p and 
niSM are the densities of the CR nuclei and gas, respec- 
tively. The diffuse synchrotron emission is proportional 
to the CR electron density, n e , and the square of the 
magnetic field strength, B, i.e., S'sync (x n P J3 2 . Assuming^ 
flux-f r eezing scaling of B oc ynisM (e.g . iRuzmaikin et al.l 
fl98J INiklas fc Bed fl997l: ICrutcherl fl999h . the syn- 
chrotron emission becomes proportional to the product 
of the electron density and the density of the interstel- 
lar gas, S'sync oc n e n\a,M- This simple scaling argument 
therefore suggests that the differences in appearance of 
the 7-ray and radio images are due to differences in the 
distributions of the CR nuclei and electrons. This was 
tested by normalizing the 7-ray and free-free corrected 
radio images by the H 1 column density map before ap- 
plying the smearing analysis, recovering results consis- 
tent with those described in 



4.2.1. Cosmic Ray Transport Properties 
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We assume that the transport of the CRs can be 
described by a simple random-walk process with an 
rigidity- dependent sp atial diffusion coefficient Dr (e.g., 
iGinzburg et al.lll980f ). The diffusion coefficient is defined 
as 



D R = D n 






(5) 



where Dq is the normalization constant and /diff is the 
characteristic distance that CRs travel after a time Tdiff. 
For our Gaussian kernels, we can relate l<as to the cor- 
responding best-fit scale lengths such that l 2 iS = a 2 = 
l 2 /2. 

If we assume that the CR electrons and nuclei are in- 
jected by the same source(s) and have been propagating 
through the ISM for the same length of time, we find that 
the ~ 20 GeV CR nuclei diffusion coefficient is ~ 4 times 
larger than that for the ~ 3 GeV CR electrons. Solving 
Equation [5] for 8, we find that the diffusion coefficient 
scales as (i?/GV) a69±015 from the exponential best-fit 
scale lengths and as (i?/GV)°' 81±0 ' 3 ° using the Gaussian 
best-fit scale lengths. Errors on 8 were estimated by a 
standard Monte Carlo approach using the uncertainties 
in the best-fit scale lengths 7 . If the minimum energy 
assumption does not hold, and we instead consider the 
extreme range of B sa 3 — 50 /zG described earlier in 
S I2.2.1[ we find corresponding 8 values of « 0.5 — 1.2. 

These model-independent values of S are consistent 
with the value of ~ 0.6 — 0.7 used in empirical diffu- 
sion models to fit the observed secondary-to-primary ra- 
tios, typically boron-to-carbon (B/C), for the Milky Way. 
In addition, our derived values for 8 generally exceed 
those for physically motivated turbulence theories, be- 
ing larger than the nominal value of | for Kolomogorov 
(|Kolmogorovl[l9 41) . and marginally consistent with the 
value of \ for Iroshnikov-Kraichnan (jlroshnikovl 119641 : 
iKraich iiau 1965|) turbulence spectra. However, the val- 
ues for the Milky Way are found for a large volume in a 
spiral galaxy, while the results of our analysis are for a 
single, highly active star-forming region: 30 Dor exhibits 
a complex network of kinematic features including slow 
(v< 100 km s" 1 ) and fast (100 km s" 1 < v < 300 km s" 1 ) 
expanding shells powered by stellar winds from youn g 
massive stars and supernovae (|Chu fc KennicuttH l994'). 

Because of strong radiative energy losses as they propa- 
gate, the corresponding diffusion lengths for the ~ 3 GeV 
CR electrons may be underestimated by our best-fit scale 
lengths, which would result in an overestimate for <5. This 
is suggested by the fact that an exponential kernel re- 
sulted in a lower minimum residual between the 1.4 GHz 
and 24 /an maps than a Gaussian kernel. However, as 
discussed below, the additional cooling of electrons may 
not significantly affect our estimate for 8. 

Taking the average derived value for 8 (i.e., 0.75), to- 
gether with the age of the current star formation ac- 
tivity responsible for supernova remnants accelerating 
the CRs, we can estimate the diffusion coefficient nor- 
malization factor Dq. The 30 Dor complex is known 

7 8 was calculated by taking 1000 random samples, picking best- 
fit scale lengths from a normal distribution having a dispersion 
set by their lower and upper-bound uncertainties. The standard 
deviation of this set is taken as the uncertainty on S. 



to contain many non-coeval stellar populations rang- 
ing in age from < 1 — 10 Myr from HST spectroscopy 
(|Walborn fc Blade s! I1997D. Rece n tly, u sing a Bayesian 
analysis. IMarti ncz-Galarza et al.l (|2011l ) fit Spitzer-IRS 
data for 30 Dor using mid-infrared SED models (contin- 
uum + lines), arriving at a luminosity- weighted age for 
the system of ~ 3 Myr. This suggests that CRs associ- 
ated with this star-forming event have been accelerated 
in the supernova remnants of very massive O-stars, and 
that the bulk of CRs have yet to be accelerated by the 
supernova remnants from the more numerous and less 
massive (i.e., ~ 8M Q ) stars with lifetimes of ~ 30 Myr. 
Note that this age is less than the estimated cooling life- 
time t coo i ~ 12 ± 2.3 Myr for the 1.4 GHz emitting CR 
electrons as they propagate through the ISM of 30 Dor 8 . 
Using this mean age of 3 Myr in Equation [5J to- 
gether with the best-fit Gaussian scale lengths corre- 
sponding to ^diff ~ 140 pc and 320 pc for the 3 GeV 
and 20 GeV electrons and nuclei, respectively, we find 
a diffusion coefficient normalization constant of Dq s=s 
0.9 — 1.0 x 10 27 cm 2 s _1 . This is more than an order 
of magnitude lower than tho se found in the Mi lky Way 
(e.g., « 5 x 10 28 cm 2 s~ 1 . iTrotta et~aT1 I20T1 . How- 
ever, analytical solutions for the perpendicular diffusion 
coefficient (i.e., diffusion across magnetic f ield lines) are 
found to be much smaller. For example, iShalchi et al.l 
(|2010f ) report a perpendicular diffusion coefficient of 
Dj_ w 3.0 - 30 x 10 26 cm 2 s" 1 . Because 30 Dor is an 
active star-forming region, which likely has a strong tur- 
bulent magnetic field and lacks a large-scale regular field 
like the Galaxy, it may not be surprising that our derived 
value for the diffusion coefficient normalization factor is 
intermediate in the range for the perpendicular diffusion 
coefficients obtained by other authors. 

4.2.2. Results Including 3-10 GeV Maps 

So far we have only discussed the results comparing the 
propagation for ~ 3 GeV CR electrons and ~ 20 GeV 
CR nuclei. While it seems that the cooling time for 
the 1.4 GHz (—3 GeV) CR electrons is a factor of — 4 
times larger than the average age of the stellar popula- 
tion in 30 Dor, trying to interpret the differences in their 
propagation lengths is complicated by the fact that the 
~ 3 GeV CR electrons lose energy more rapidly than the 
~ 20 GeV CR nuclei. That the electron energy losses are 
important is suggested by the result that the residuals be- 
tween the 1.4 GHz and 24 /j,m maps are decreased more 
when using exponential kernels rather than Gaussian ker- 
nels. We therefore include the 3—10 GeV 7-ray map to 
provide an additional estimate for the propagation length 
of CR nuclei because it is sensitive to much higher en- 
ergy CR nuclei (i.e., ~ 70 GeV for a proton energy index 
of 2.1). Employing these data, we can try to determine 

8 This maximum lifetime is estimated using the combined energy 
losses due to synchrotron, inverse Compton, bremsstrahlung, and 
ionization processes (see, e.g., Murphy 2009). Using the previously 
estimated values for the minimum energy magnetic field strength 
(B m ; n ss 11 ± 4/iG; see §2.2), the radiation field energy density 
(Umud ~ 2.39 ± 0.17 X 10 — 12 erg cm -3 ; see §2.4), and assuming an 
average ISM density of nigM ~ 2±1 cm~ 3 (Kim et al. 2003), we es- 
timate a cooling time for the observed 1.4 GHz emitting (~ 3 GeV) 
electrons of ~12 ± 2.3 Myr. We note that the individual lifetimes 
estimated against synchrotron, inverse Compton, bremsstrahlung, 
and ionization losses are 33, 55, 43, and 140 Myr, respectively 
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whether the difference between the CR electrons and pro- 
tons is in fact the result of rigidity-dependent diffusion. 

Applying the same smearing analysis to the 1.4 GHz 
radio, 1 — 3 GeV, and 3—10 GeV 7-ray maps, but limiting 
the area of the residual calculation to that where there 
are 3 — 10 GeV 7-ray events (i.e., w \ the area having 
1 — 3 GeV 7-ray events), we obtain the following results. 
The best-fit exponential and Gaussian scale lengths are 
50 and 150 pc, respectively, for ~ 3 GeV electrons, 150 
and 300 pc, respectively, for ~ 20 GeV CR protons, and 
200 and 400 pc, respectively, for ~ 70 GeV protons. Be- 
cause the detected 3—10 GeV events are concentrated 
around the center of 30 Dor, lacking a significant ex- 
tended component as observed in the 1 — 3 7-ray map, it 
is not surprising that the scale lengths obtained are gen- 
erally smaller than those calculated from the residuals 
over a larger area. 

Using all six possible combinations of these values 
to solve for 5 results in a median and dispersion of 
5 « 0.66 ± 0.23. This is consistent with the energy de- 
pendence estimated above using the 1.4 GHz radio and 
1 — 3 GeV 7-ray maps. We note that the value of 5 esti- 
mated by comparing only the 7-ray maps, which is inde- 
pendent of assumptions for the magnetic field strength 
and proton energy index, is ~ 0.51, whereas the me- 
dian when only including the 1.4 GHz maps is w 0.79, 
suggesting that the best-fit scale lengths for the CR elec- 
trons may be underestimated due to additional energy 
losses. This is consistent with exponential kernels work- 
ing slightly better to tighten the correlation between the 
1.4 GHz and 24 /mi maps compared to Gaussian kernels. 
Additionally, the estimate for the value of 5 using only 
the 7-ray maps is independent of assumptions for the 
magnetic field strength However, we emphasize that the 
3 — 10 GeV 7-ray map is statistically limited, having > 3 
times fewer events than the 1 — 3 GeV maps. 

4.2.3. Propagation Length vs. Star Formation Activity 

Having estimated the average distance traveled by CR 
electrons for an individual star-forming region, 30 Dor, 
it is interesting to see how these results compare with 
similar estimates of CR electron propagation distances 
for entire galaxies. iMurphv et "all (|2006bL 120011 reported 
a correlation between the CR electron propagation dis- 
tance and galaxy surface brightness such that CR elec- 
trons are found to travel shorter distances, on average, 
in galaxies having higher star formation a ctivity. 

In Figure [5] we plot the results from IMurphv et al.l 
(|2008l i.e., the best- fit exponential scale lengths pro- 
jected in the plane of the sky versus t/ ra d) along with our 
estimates for 30 Dor. We note that various systematics 
and selection effects (e. g., distance, inc l inatio n) were rig- 
orously investigated by M urphy et al.l (|2008[) . They ad- 
ditionally tested the differences between isotropic kernels 
and kernels projected in the plane of each galaxy disk, 
finding a preference, i.e., smaller residuals, for isotropic 
kernels. 

The C/ ra d estimate for 30 Dor was calculated from the 
surface brightness measured over an area of 2.4 kpc 2 , 
the same region for which the residuals were estimated 
(see fc|2.4[) . For this case, the best- fit scale length was 
measured by comparing the free- free corrected 1.4 GHz 
with smoothed 70 /x m maps to allow a pro per comparison 
with the results of IMurphv et al.l ([2008). The best-fit 
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Figure 5. Best-fit exponential scale lengths plotted against radi- 
ation field energy density (measu red by the infrared su rface bright- 
ness) for the galaxies studied in Murphy et al. (2008, diamonds), 
the LMC (asterisk), and 30 Dor (filled star). The open diamonds 
indicate spirals, while the filled dia monds show the loca tion of 
dwarf irregular galaxies included in the Murphy et al. (2008) study. 
The location of 30 Dor (filled star) follows the trend of decreasing 
exponential scale length, indicating shorter propagation lengths of 
CR electrons with increasing radiation field energy density for the 
more active star-forming disk galaxies in the Murphv et al. (2008) 
study. We plot the CR electr on propagation distance for 30 Dor 
as an upper limit (see i|4.2,3l ) For the dwarf irregular galaxies, 
the short scale lengths are thought to be due to an increase in the 
escape of the CR electrons from the systems. If we set the LMC 
at the mean distance of the Murphy et al. (2008) dwarf irregulars 
(i.e., ~ 3.9 Mpc), its propaga tion scale length is more consistent 
with the Murphv et al. (2008) spirals. 
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Figure 6. The same as Figure \E\ except that we show only 
galaxies having high disk-averaged star formation rates. A distinct 
trend of decreasing CR electron propagation length with increas- 
ing radiation field energy density can be seen. Galaxies from the 
Murphv et al. (2008) study have been labeled. The dotted line is 
an ord inary least squares fit only to the galaxies from Murphy et al. 
(2008). The locations of 30 Dor (filled star, as an upper limit) and 
the entire LMC (asterisk; scaled to d = 3.9 Mpc) are plotted. 

exponential and Gaussian ker nel scale length s were and 
50 pc, respectively. Because IMurphv et al.l (|2008l ) plot 
exponential kernel scale lengths, we plot the result from 
smoothing the 70 fjjn map of 30 Dor with a Gaussian 
kernel as an upper limit. 

Included in Figure [5] are both spirals and star-forming 
irregulars. The latter have markedly small best-fit scale 
lengths for their values of C/ rac j relative to the sample 
of spirals. This is thought to be caused b y increased es- 
cape of CR ele ctrons in these syst ems (e.g.. lCannon et ahl 
120051 120061: IMurphv et all 120081 ). To see how the be- 
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havior of entire galaxies compares with the LMC as a 
whole, we repeat the image smearing analysis using the 
70 /jm maps aft er first project i ng the LMC to the mean 
distance of the iMurphv et al.l (|2008[ ) sample irregulars, 
d = 3.9 Mpc to take possible resolution effects into ac- 
count. We calculate the residuals and C/ rac i within an 
aperture having a radius of « 3.1 kpc (2.'75 at the dis- 
tance of 3.9 Mpc) centered at 05 h 18 m 50 s , -68°50'49" 
(J2000). This is approximately the area over which the 
24 /im emission was detected at > 3cx level. The best- 
fit exponential and Gaussian scale lengths are 200 and 
400 pc, respectively. This places the LMC close to the 
main trend among the star-forming spirals. 

Focusing only on the region containing the star- 
forming spirals, we re-plot the disk galaxies along with 
30 Dor and the LMC in Figure El The do tted-line is a 
least-squares fit to the lMurphv et all (J2008I ) galaxy disks, 
excluding the LMC and 30 Dor. Clearly, the position 
of 30 Dor in the plot is consistent with the empirical 
trend describing spiral galaxies, but extrapolated to a 
surface brightness greater by almost an order of magni- 
tude. Thus, this scaling relation appears to operate on 
the scales of entire galaxies all the way down to individ- 
ual star-forming regions. 

5. CONCLUSIONS 

Using a phenomenological image smearing model, we 
have estimated the typical propagation length of ~ 
3 GeV CR electrons and ~ 20 GeV CR nuclei by com- 
paring the spatial distributions of smoothed 24 /im and 
(free- free corrected) 1.4 GHz and 1 — 3 GeV 7-ray maps, 
respectively. Below we list the major results and conclu- 
sions: 



value is consistent with that obtained by including 
the more statistically limited 3 — 10 GeV 7-ray map 
(~ 70 GeV CR protons; S « 0.66 ± 0.23). 

4. The value of S reported here is larger than 
the spectral index of the diffusion coefficient 
assuming Kolomogorov turbulence and Milky 
Way secondary-to-primary ratios, and marginally 
consistent with Iroshnikov-Kraichnan turbulence. 
This may reflect the fact that 30 Dor region ex- 
hibits complex kinematic features, and fast expand- 
ing shells, resulting in a larger value of 5. 

5. Assuming that the CRs in 30 Dor are as old as 
the average stellar population (~3 Myr), we esti- 
mate a diffusion coefficient normalization constant 
of D « 0.9 - 1.0 x 10 27 cm 2 s- 1 . This value is 
less than and similar to mo del- dependent estimates 
of the parallel and perpendicular diffusion coeffi- 
cient for the Galaxy, respectively. The similarity 
between our estimate and perpendicular diffusion 
coefficient values may be due the 30 Dor region 
having a magnetic field configuration that is highly 
turbulent. 

6. We place our results for 30 Dor, and the LMC as 
a whole, in the context of the scaling relation be- 
tween the typical CR electron propagation length 
and disk-averaged star formation activity for en- 
tire disk galaxies, where CR electron propagation 
is found to decrease with increasing star formation 
activity. This relation appears to apply to the LMC 
and the individual star-forming region of 30 Dor. 



1. We estimate the typical distances traveled by 
-3 GeV CR electrons from 30 Dor to be - 100 - 
140 pc. This is factor of ^2 smaller than that es- 
timated for ~20 GeV CR nuclei, which is found to 
be - 200 - 320 pc. 

2. In our image-smearing analysis, we find that ex- 
ponential kernels work slightly better to tighten 
the correlation between the 1.4 GHz and 24 /xm 
maps compared to Gaussian kernels. In contrast, 
both exponential and Gaussian kernels are found 
to work equally well to tighten to correlation be- 
tween the 1 — 3 GeV 7-ray and 24 /xm maps. This 
difference suggests that, unlike the CR nuclei, CR 
leptons suffer additional energy losses as they prop- 
agate through the ISM near 30 Dor on timescales 
less than, or comparable to, the diffusion timescale. 

3. Assuming that the CR electrons and nuclei are pro- 
duced by the same sources, and that their prop- 
agation is well described by a random walk, the 
differences in their estimated propagation lengths 
suggest differences in their associated spatial diffu- 
sion coefficients. This allows us to make the first 
model-independent measurement of the energy de- 
pendence of the diffusion coefficient for an external 
galaxy. For CRs produced in this star-forming re- 
gion, the ~20 GeV CR proton diffusion coefficient 
is ~4 times larger than that for ~3 GeV CR elec- 
trons, scaling as (R/GV) S where d w 0.7-0.8. This 
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